Electromagnetic Energy, Momentum, and Angular Momentum in an Inhomogeneous 

Linear Dielectric 



o 

(N 



(N 



Michael E. Crenshaw and Thomas B. Bahder 

US Army Aviation and Missile Research, Development, 

and Engineering Center, Redstone Arsenal, AL 35898, USA 

(Dated: March 1, 2013) 

In a previous work, Optics Communications 284 (2011) 2460-2465, we considered a dielectric 
medium with an anti-reflection coating and a spatially uniform index of refraction illuminated at 
normal incidence by a quasimonochromatic field. Using the continuity equations for the electro- 
magnetic energy density and the Gordon momentum density, we constructed a traceless, symmetric 
energy-momentum tensor for the closed system. In this work, we relax the condition of a uniform 
index of refraction and consider a dielectric medium with a spatially varying index of refraction 
that is independent of time, which essentially represents a mechanically rigid dielectric medium due 
to external constraints. Using continuity equations for energy density and for Gordon momentum 
density, we construct a symmetric energy-momentum matrix, whose four- divergence is equal to a 
generalized Helmholtz force density four-vector. Assuming that the energy-momentum matrix has 
tensor transformation properties under a symmetry group of space-time coordinate transformations, 
we derive the global conservation laws for the total energy, momentum, and angular momentum. 
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I. INTRODUCTION 

Starting with the hydrodynamic continuity equation, 
Umov [1] obtained an expression for energy continuity in 
a continuous spatial flow of an electromagnetic field in 
1874 0-01 • A decade later, Poynting [2[ derived a sim- 
ilar energy continuity equation as a general theorem of 
the macroscopic Maxwell equations. Poynting's theorem 
is generally preferred because it can be derived directly 
from the macroscopic Maxwell equations 
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of classical continuum electrodynamics. Here, the 
Maxwell equations are written in Heaviside-Lorentz units 
for a nonmagnetic linear medium in the absence of free 
charges and currents. We have assumed that the index 
of refraction, n = n(r), depends on position, occupies a 
finite region of 3-dimensional space, and is independent 
of time. Poynting's theorem can be derived by subtract- 
ing the scalar product of the Maxwell- Ampere law, Eq. 
(|1.1|) . with E from the scalar product of Faraday's law, 
Eq. (fl~2j) . with B to obtain 

n 2 dE „ 10B „ 
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Upon application of a vector identity and the definition 
of the energy density, 



Pe = i(n 2 E 2 + B 2 ), 



(1.6) 



the preceding equation becomes Poynting's theorem 

V-c(ExB) = (1.7) 



dt 



and defines Poynting's energy flux vector Sp = c(E x B). 
In a previous work [5|], we considered a dielectric 
medium with an index of refraction that was time- 
independent, spatially uniform, and covered with a thin 
gradient-index antireflection coating. For this case, we 
showed that the energy 
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and the Gordon [6j momentum 

G G = ~(riE x B)dv 
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(V x E) B- (V x B) -E = 0. 
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are the conserved electromagnetic quantities with inte- 
gration performed over all space a. We constructed the 
corresponding traceless, symmetric energy-momentum 
tensor, whose four-divergence provided the continuity 
equations for the energy and momentum densities of the 
closed system. This result incorporated a condition of 
Vn/n <C 1/A corresponding to unimpeded flow in the 
absence of external forces. 

In the current work, we consider the case of an in- 
homogeneous dielectric medium in which the condition 
S/n/n <C 1/A no longer holds. The assumption that the 
index of refraction is independent of time necessarily im- 
plies that the dielectric medium is mechanically rigid Q 
and that no momentum is transferred to the dielectric 



medium. In other words, at each spatial point, there 
is effectively an external field that acts as a constraint 
that holds the dielectric in place. Therefore, the sys- 
tem that we are considering is not a closed system and 
the continuity equations and global conservation equa- 
tions will reflect this feature. The external field that acts 
as a constraint is not a field that we impose, instead it 
arises as a result of the way in which we arrange terms in 
the continuity equations. We construct the energy and 
momentum continuity equations and find that the spa- 
tial gradient of the refractive index appears in the con- 
tinuity equations in a generalized Helmholtz force den- 
sity. This generalized Helmholtz force density is a spa- 
tial and time- varying field that provides the constraint 
for the dielectric medium to be mechanically rigid. We 
then write the continuity equation as a four divergence of 
a symmetric energy-momentum matrix. Assuming that 
the energy-momentum matrix has tensor transformation 
properties under a symmetry group of space-time coordi- 
nate transformations, we apply the four-dimensional di- 
vergence theorem to derive global conservation laws for 
the total energy, momentum, and angular momentum. 
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The Gordon momentum in Eq. (|2.2[) has been scaled by 
c so that the equations are in the same units. Next, we 
apply the vector identity 



V x (tpa) = Vt/> x a + tpV x a 



(2.3) 



to Faraday's law and write the macroscopic Maxwell 
equations, Eqs. (11.11) and (11.21) . as 
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The variant form of Maxwell's equations, Eqs. (|2.4|) and 
(I2.5[) . are mathematically equivalent to the original ver- 
sions, Eqs. (jl.ll) and (|1.2p . respectively. Substituting the 
Maxwell equations, Eqs. (f!Q|) and (|2T5j) . into Eqs. (|2~Tj) 
and (|2.2|) . we produce the energy continuity equation 



II. ENERGY AND MOMENTUM CONTINUITY 
EQUATIONS 



!^ +V -( n ExB) = --(nExB) (2.. 

c at n 

and the (Gordon) momentum continuity equation 



The energy-momentum tensor is a concise way to rep- 
resent the local continuity of energy and momentum 
of a field. While simple in concept, the form of the 
energy-momentum tensor for the electromagnetic field 
in a dielectric has been at the center of the century- 
long Abraham-Minkowski controversy @. The tensor 
that was proposed by Minkowski [9j in 1908 is not diag- 
onally symmetric, a fact that is adverse to conservation 
of angular momentum |10| . Abraham 11 1 subsequently 
proposed a symmetric tensor at the expense of a phe- 
nomenological force. The disagreeable properties of the 
Minkowski and Abraham energy-momentum tensors are 
manifestations of underlying conservation issues: neither 
the Minkowski momentum nor the Abraham momentum 
is conserved. In this section, we derive the energy and 
momentum continuity equations for an inhomogeneous 
dielectric and construct the corresponding tensor con- 
tinuity equation and energy-momentum tensor. In the 
following section, we will obtain the globally conserved 
quantities — the total energy, the total momentum, and 
the total angular momentum system, from the energy 
momentum tensor. 

Starting from the Maxwell Eqs. (|l.l| )- (|1.4|) . we use the 
continuity of the electromagnetic energy density and the 
Gordon momentum to obtain a symmetric stress-energy 
tensor. We begin with the temporal derivatives of the 
energy density and the Gordon momentum density, 
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Here, we have used E q. Ill -3D and the definition of the 
Maxwell stress tensor [lfj, [l2| 



Wij = [ -nE.nEj - B.B, + - (nE • nE + B • B) % 
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We can write the energy continuity equation, Eq. 
and momentum continuity equation, Eq. (J2.7I) . as the 
matrix differential equation 
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with summation over repeated indices. The quantities 
that appear in Eq. (12.91) are a four-divergence operator 
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the array 



rpafi 
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and a generalized force density four-vector 
r = (Vn-(ExB),-fo), 
where 



(2.12) 



fe = -nE 2 Vn + 2nE (E x Vn) + n 2 E(V • E). (2.13) 

Using Eq. (11.31) . the last two terms on the right side of 
Eq. (|2.13p cancel, so that 



f ff = -nE 2 Vn 



(2.14) 



in the absence of free charges. 

The array in Eq. (|2.1ip appears to have the properties 
of an energy-momentum tensor. By construction, the 
operator defined in Eq. (J2.10I) applied to the rows of the 
array in Eq. (|2.11[) generates continuity equations for the 
electromagnetic energy and the momentum. The same 
operation applied to the columns 

B a T afj = f, (2.15) 

generates the same continuity equations by symmetry 

T a(3 _ T (3a_ ^jg) 



The array has a vanishing trace 



(2.17) 



corresponding to massless particles [10|, [l2[ . 

The new feature of this result is the appearance of 
the four- vector f a in the continuity equation (|2.9I) . The 
appearance of the four- vector f a in the divergence of the 
stress energy tensor in Eq. (I2.9P is a result of the fact 
that the system is not closed. When the gradient of the 
index of refraction is non-zero, there is a back-action on 
the field altering its spatial properties, so the field does 
not experience "unimpeded flow" . It is possible to define 
an effective stress-energy tensor that takes into account 
the back-action on the field by introducing 
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The continuity of energy and momentum can then be 
expressed by 



a/3 , ,a{}\ _ 



d a (T ap + t ap )=Q 
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where the tensor t a ^ contains the influence of the inho- 
mogeneity of the material and is zero for homogeneous 
dielectrics. 



III. GLOBAL CONSERVATION EQUATIONS 

If we assume that T a/3 transforms as a tensor under 
some symmetry group of space-time coordinate transfor- 
mations, then we can apply the four-dimensional diver- 
gence theorem to obtain global conservation equations. 



Integrating Eq. (J2.9I) over a four- volume, dH, = cdt d 3 x = 
cdt dv (where i — t/ri), between hypersurfaces [Ty] of 
constant time at i\ = ti/n and t% = tz/n, we have 



*2 rt 2 

qO 



dpT aV cdt dv = f a cdt dv, (3.1) 
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where the integrals over dv — d 3 x are 3-dimensional vol- 
ume integrals over the volume containing the field. Ap- 
plying the four-divergence theorem results in 



(T a0 (t 2 ) - T Q °(ti)) dv = c / -f a dt dv. (3.2) 

n 



When a = 1, 2, 3, Eq. fl3J2]) reduces to 

G(i 2 ) - G(ti) = I E 2 Vn dt dv (3.3) 

where G(t) is the Gordon momentum at time t, and the 
time integration is between t = t\ and t = t%. Equa- 
tion (J3.3I) shows that the Gordon momentum of the field 
is not constant, i.e., there is a source or sink of momen- 
tum provided by the external constraint field f#. For the 
case a = 0, Eq. (|3.2[) gives 

f (Pe(t 2 ) - Pe{h)) dv= J -^■(ExB)dtdv (3.4) 

where p e is given by Eq. (|L .€>[) . The left side of Eq. p.4p 
is the difference in total energy at two different times, 
U{ti) — U(t\). We see that a non-zero spatial gradient in 
the index means that the field can gain or loose energy. 
The results given in Eqs. (J3.3I) and (|3.4p assume that the 
matrix T a/3 in Eq. (|2 .9[) transforms as a tensor under 
some symmetry group of space-time coordinate transfor- 
mations. Note that when Vn = 0, Eqs. p. 31) and (|3.4p 
show that the Gordon momentum G(i) and the energy 
U(t) are constants. 

In recent years, there is intense interest in angular mo- 
mentum carried by the electromagnetic field, see the re- 
cent review and references cited therein [lj, [lj|. How- 
ever, the angular momentum carried by the electromag- 
netic field in a dielectric environment is no less unsettled 
than the linear momentum case [Tfj. We can define the 
four-tensor of angular momentum density in terms of our 
energy-momentum tensor as jlfjj] 
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Continuity of angular momentum is given by 
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The divergence of m Q ^ 7 is not zero, thereby indicating 
that there is a source or sink of angular momentum den- 
sity, due to the gradient of the index of refraction. 

Once again, we assume that there exists a symmetry 
group of space-time coordinate transformations, so that 



m Q/37 is a tensor. As above, we can then use the four- 
divergence theorem to obtain 
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(<i)) dv 



1 



(x a f-x f a )dtdv 

where the time integral is between the two times t\ and 
t 2 and where the volume integral dv is over the portion of 
three-dimensional space containing the field and includes 
the region where n(r) > 1. When a and /3 take values 
i,j = 1, 2, 3, we have 



^i Q t+\ — rr-ifJ 



(t)=x l P-xif 



(3.8) 



where p are the components defined in Eq. (|2.14p . Equa- 
tion (J3.7I) gives the change in total angular momentum 
of the field, AM, between time t\ and t<2, and can be 
written as 

AM = / (m(t 2 ) - m(ti)) dv = I E 2 (r x Vn) dtdv 

where the m(£) is the angular momentum density at po- 
sition r at time t, defined by 

m(t) = r x g G (t) (3.10) 

and is defined in terms of the Gordon momentum density 



gG (i)=nExB/c. 



(3.11) 



Equation (|3. 91) shows that when r x Vn is non-zero, then 
the total angular momentum of the field (as expressed 
through the angular momentum density) can change. 
Once again, we remind the reader that we have assumed 
the index of refraction as isotropic, constant in time, but 
varying in position. Equation (|3.9[) essentially shows that 
a particular spatial distribution of r x Vn can lead to a 
back-action on the field that can alter the field angular 
momentum. Indeed this has been exploited in a number 
of recent experiments [13] . 



IV. SUMMARY 



In a previous work, we considered a dielectric medium 
with an index of refraction that was time-independent 
and spatially uniform, and the dielectric was covered with 
a thin gradient-index antireflection coating Q . We found 
that the total energy and the Gordon momentum were 
conserved quantities (constant in time). In this work 
we have relaxed these conditions to include a dielectric 
medium that has a spatially varying index of refraction 
that is constant in time. The fact that the index of re- 
fraction is constant in time essentially means that the 
dielectric medium is mechanically rigid and subject to 
an external constraint. This constraint means that we 
are not dealing with a closed system. Using the continu- 
ity equations for energy density and momentum density, 
we derived a symmetric energy-momentum tensor for the 
electromagnetic field. Due to the fact that the system is 
not closed, the divergence of the stress-energy tensor is 
not zero, but equal to a generalized Helmholtz force vec- 
tor that represents the constraint. Similarly, the diver- 
gence of the angular momentum density is not zero, due 
to the external constraint of a time independent index of 
refraction. We found that the total energy, Gordon mo- 
mentum, and total angular momentum are not conserved 
because of the external constraint on time independence 
of the index of refraction. However, the time dependence 
of the total energy, Gordon momentum, and total an- 
gular momentum is related to the constraint, which is 
proportional to the gradient of the index of refraction. 

Finally, we note that time t has been renormalized 
to i = t/n(r) in the intermediate steps used to obtain 
the symmetric energy-momentum tensor and its conti- 
nuity relation given by its four-divergence in Eq. (12.91) . 
However, the final conservation laws of total energy, to- 
tal momentum and total angular momentum, given by 
Eqs. (|3.4[) . (|3.3p . and (|3.9[) . respectively, are expressed in 
terms of the unrenormalizcd time t. 
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